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Abstract 
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Recently some authors have broadened the scope of canonical quantum mechanics 
by replacing the conventional Hermiticity condition on the Hamiltonian by a weaker 
requirement through the introduction of the notion of pseudo-Hermiticity. In the 
çSj . present study we investigate eigenvalues, transmission and reflection from complex 

optical potentials enjoying the property of pseudo-Hermiticity. 
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Bender, his collaborators and others [1-23] in a series of papers, investigated 
some non-Hermitian Hamiltonians which violate parity (P) and time reversal 
(T) symmetry but are PT invariant. These systems comprised of a particle 
moving in a complex potential, the real part of which is parity-even while 
the imaginary part is odd. Through various examples it was found that the 
energy eigenvalues were real and bounded from below. Though in some cases, 
for a range of parameters (contained in the potential), complex eigenvalues do 
occur, but these are associated with a spontaneous breaking of PT symmetry. 
Subsequently Mostafazadeh [24] has provided a bàsic mathematical setting for 
such systems by introducing the notion of pseudo-Hermiticity of the Hamil- 
tonian H via the condition W = rjHr]^ 1 where i] is a Hermitian, linear and 
invertible operator. It is suitable to define an indefmite inner product of two 
state vectors: <^ifji\ipi2^> v = (^íl^l^) which is time translationally invariant 
with rj playing the role of a mètric operator in Hilbert space. One can then go 
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on to show that the eigenvalues of H are either real or occur as complex con- 
jugate pairs and the eigenvectors constitute a complete biorthonormal system. 
For the examples considered by Bender and others rj is the parity operator 
('I /')• 

In the present study we attempt to put such Hamiltonians in a more 'physical' 
setting and moreover extend the discussion to continuum states as well. We 
encounter some rather amusing features. We recali [25] the optical (or cloudy 
crystal ball) model of the atòmic nucleus, where the interaction of a neutron 
(or some other projectile) with the nucleus is described through a complex 
potential V + iW, where V and W are its real and imaginary parts respec- 
tively. From the continuity equation obeyed by the probability and probability 
current densities (following from Schròdinger equation) it is easily seen that 
non-vanishing W implies a sink (or source) for the probability, depending on 
the sign of W, which is taken to correspond to absorption (or emission) of 
partides with respect to the incident beam, and this furnishes a rather use- 
ful phenomenological description of elàstic scattering in the presence of open 
inelastic channels. In one of the versions of this picture, known as the sur- 
face absorption model, the imaginary part (W) of the potential is taken to 
be proportional to the spatial derivative (gradient) of the real part (V). We 
shall see that a one-dimensional 'cartoon' of the surface absorption version of 
the optical model leads naturally to a PT- invariant (or more appropriately a 
pseudo-Hermitian) Hamiltonian. We illustrate some novel features of pseudo- 
Hermitian Hamiltonians through two examples so chosen such that the results 
are expressible in terms of elementary functions. 

Consider a particle of mass m moving in one dimension under the influence of 
a potential the real part (V) of which is an attractive square well of depth Vq 
and range a and whose imaginary part (W) is proportional to the derivative 
of the real part. Accordingly the motion is governed by the Hamiltonian 

H = L + v + m = L - V K ~ «> + < mx ~ i» " s(x + > (1) 

Note that the Hamiltonian is not parity invariant as V + iW — >V — iW 
(W being an odd function of x ), nor is it time reversal invariant as the 
operation of complex conjugation takes V+iW — >V—iW. However the system 
enjoys symmetry under the combined PT transformation. More relevantly the 
Hamiltonian is pseudo-Hermitian viz. W = PHP^ 1 with the parity operator 
P being allowed to play the role of the 'mètric' operator (77) in the Hilbert 
space. 

Analysing the bound state situation to begin with, it is convenient to put 
the energy eigenvalue E = —B with B > and to define = j3 2 , so that 
ipi(x) = Ae^ x for x < — | (as ip must vanish as x — >— 00) and ipjjj(x) = De~P x 
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in the region x > +| (as ip must vanish as x — > + oo). In the range \x\ < | we 
have the solution íjjn(x) = Cie tqx +C2e~ igx where q 2 = 2m ^~ B ^> , Implementing 
the continuity and jump conditions i^i(^) = Í>ii(^y) = ^iii(^), 

fn(=r) - Vi(=?) = -àM=r), fmi^) ~ ^ni^f) = iMni^) where 
À = we arrive at the eigenvalue condition 

(q 2 -[3 2 -\ 2 ) Sin(qa) = 2(3q Cos{qa). (2) 



Before going on to a discussion of the solution to this equation, it is revealing 
to ask the question: Is there a minimum depth Vq below which a bound state 
(with the usual meaning of the concept as elaborated below) does not exist ? 
For the square well potential (with À = 0) we know that there is at least one 
bound state howsoever weak the attraction ( given by Vo ) may be. Also for 
^j 2mVoa 2 < | ^here j s on iy one bound state. Let us choose the real part of the 

potential to lie in this region ( say 2m ^ a2 = 1 ) and look for the effect of the 
imaginary part (À) on the binding energy by numerically solving eq.(2). The 
result is shown graphically in Fig.l(a). Observe that the binding energy goes 
to zero at À = = 1- This is readily seen by looking at the zero binding 
condition, namely, (3 = in the eigenvalue equation [eq. (2)] which yields the 
relation 



2mÀ 2 .^. , 2mV . „ 
(V - -^) S ^-^ a ) = 0- (3) 



The relevant result for this minimum depth is thus 

T, 2mX2 í A\ 

Vo = (4) 



In this limiting situation the wavefunction corresponding to zero binding (B = 
0) is 

ipi = A for x < -§, 

■ mXa ■ 2mXx 

$ n = Ae' 1 ^- 1 -^- for -§ < x < +§ 

■ 2m\a 

i\)ui = Ae for x > +§. 

To examine analytically what happens when the strength of imaginary part 
crosses the critical value, let us choose for simplicity \j 2m ^ a ' 1 = | and look 

in the neighbourhood of the point À 2 = 2r ?¥° — ^ 2 vi z ^ 2fl2 — oj 2 a 2 — e 
where, with e positive, a bound state should exist. With sufficiently small e 
the roots of the transcendental equation [eq. (2)] reduces to finding the zeros 
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of a cúbic form and leads us to three solutions (3a = — |e ± + 0(e 2 ) 

and /9a = — 2 + O(e). In order to have ip(x) = e~^ x — > as ï — > oo 
and ip(x) = e@ x — > as x — ► — oo corresponding to bound states it is 
necessary to have (3 > and thus only one of these three roots is physical 
namely (3a = — |e + + 0(e 2 ) , the other two being unphysical. With À 
vàlues beyond the critical value looking again at the neighbourhood of that 

point we obtain /3a = §|e| ± i\f^ and the third root at —2 + 0(e) which as 
before is inadmissible. While the wavefunctions corresponding to the complex 
conjugate roots are well behaved in the sense that |-0(x)| 2 — ► as x — ► ±oo 
and are square integrable, nevertheless there is a serious difficulty regarding the 
physical interpretation of the time dependence of the probablity. Thus while 
the wavefunction corresponding to the root with negative imaginary part can 
be thought of as a decaying state, its inevitable partner with positive imaginary 
part grows with time and is thus physically unacceptable. Moreover these two 
states are not orthogonal in the sense that <C ipi\ip2 ^»7= (V'iMV^) 7^ 
(where rj — P ) while <C ^ v = ip2\ "02 Thus the system itself 

ceases to have usual physical significance when these complex conjugate roots 
make their appearence, namely, when the strength of the imaginary part of the 
potential À exceeds the critical value. It may be noted more generally that if 
the real part of the potential were deeper so as to support more than one bound 
state then the boundstates continue to be real (and interpretable) as long as 
the imaginary part is less than the critical value at which the least bound 
state complexifies and the system no longer sustains its physical meaning. 
It may also be noted that in the region of parameter space when complex 

energy eigenvalues occur, the wavefunctions are given by ip(x) ~ e 8 a V 2 a 

_3| e |£_j. /Hi 

for x < and i((x) ~ e 8 a V 2 a for x > (for À slightly greater than 

2 V -\e\-—i\h-^-- 

the critical value -jp-) [as also the other root giving tp{x) ~ e 8 a V 2 a for 

3 I c I X 1 £ / I e I x 

x < and ip(x) ~ e 8 a V 2 a for x > 0]. These are not FT-invariant and 
as such with the Hamiltonian enjoying the symmetry and the states violating 
it, we see that PT is spontaneously broken. 



This situation may be contrasted with the system of a particle moving in a 
complex Morse potential in one dimension which was so contrived as to have 
only real energy eigenvalues [26]. It has also been pointed out [27] that for 
this potential the resulting Hamiltonian is ?]-pseudo-Hermitian with rj = e~ 6p 
where p is the momentum operator. As has been shown by Mostafazadeh [28] 
the necessary and sufficient condition that a pseudo-Hermitian Hamiltonian 
has only real eigenvalues is that we may write 77 = A^A, where A is a linear 
and invertible operator. Clearly for the operator e~ ep the corresponding A is 
e ~2 e P. Noting that this is not the case with the examples chosen here, we have 
real eigenvalues only for a certain regime of parameters. 
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A convenient framework for the discussion of transmission and reflection for 
the problem at hand is provided by the S— matrix approach. We introduce 
asymptotic channel states \m,k), with m = R, L : (x\R,k) = e lkx and {x\L, k) = 
e~ lfcr , where R and L stand for right moving and left moving free particle states 
and the wave-number k = ^Jj^E. One can then define the S operator whose 
on-shell matrix elements (m,k\S\n,k) = S mjTl (k) gives the probability am- 
plitude for a state starting off in the remote past as \n,k), to be found, as 
a result of evolution through the interaction with the potential, in the state 
\m, k) in the remote future. This 2x2 matrix S would have been unitary if the 
potential were real, but in the present case this will not be so. With the use 
of these conventions for enumerating channels, it is evident that the S matrix 
elements are related to the familiar transmission (t R and ti) and reflection (r R 
and tl) amplitudes for right and left travelling partides (indicated through 
the subscripts): 



(Srr S rl \ _ ít R r L \ 

\S L r S ll ) ~ \r R t L ) [ ' 

With hermitian Hamiltonians the states evolve in a unitary manner and the 
S— matrix obeys the unitarity condition S^S = I which in this case would 
imply the relation \t R \ 2 + |r fí | 2 = 1, |íl| 2 + \r L \ 2 = 1 and t* R r L + r R t L = 0. The 
first two of these conditions are nothing but the conservation of probability 
for the right and left incident beams respectively while the third describes the 
phase relationships. Note that in a left-right symmetric situation this reduces 
to t*r + r*t = or that the transmission and reflection amplitudes are out of 
phase by |. In our case since the Hamiltonian is P-pseudo hermitian, it is clear 
that as W = PHP^ 1 the S— matrix obeys the pseudo unitarity condition 
P^^^S^PS = I. The P operation in the \L) , \R) basis in which we have 
expressed the S— matrix is given by the matrix 



(6) 



This implies that t* L t R + r* L r R = 1, r* L t L + t* L r L = and r* R t R + r R t* R = 0. The 
last two conditions imply that the reflection and transmission amplitudes are 
out of phase by | for both the left and right incident beams. 

For the particular model defined by the Hamiltonian given in eq.(l), the trans- 
mission amplitudes for left and right incident beams are explicitly found to be 
the same, and in our example is given by 

t =t = 2 V ke ~ lka m 

R L 2qk Cos(qa) - i (ç 2 + k 2 - À 2 ) Sin(qa) ' 
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The fact that Ír = Íl follows from PT symmetry of the potential. Since the 
Hamiltonian is PT— symmetric, hence the S— matrix is also PT— symmetric. 
For a hermitian Hamiltonian which is invariant under the time-reversal it is 
well known (see for example the discussion in the text-book by A.S. Davydov 
[29]) that the S— matrix gets transposed under the operation of time reversal 
(T). Following the same procedure it can be seen that even though, as in 
our case, the Hamiltonian changes under time reversal (T) and parity (P) 
separately , but is symmetric under their joint operation (that is PT) the 
S— matrix suffers a transposition under T. Accordingly (PT)H(PT)^ 1 = H 
while PHP- 1 = i/t ; and thus THT^ 1 = PHP- 1 = H\ Introducing a unitary 
operator O (0^0 = i) such that OH* = WO then following Davydov [29] wc 
can assert that T = OK, where K is the operator for complex conjugation. 
Thus the wavefunction corresponding to the time reversed state is tp- a = 
Tip a = OKipa = Oip* a . The S— matrix element between the time reversed 
states is 

S- a ,- b = (1>- a \S\il>- b ) = (T^ a \S\T^ b ) = (Or a \S\Oip b *) 

= moiso\rt) = (r a \(&y\iit) = s ba 

where we have used the relation OH* = H^O to obtain O^SO = (S't)* an d 
then the fact that (^|(5t)*|^> = (Al&lA)* = (4>b\S\4> a ) = S ba . 

Therefore the S— matrix get transposed under the operation of time rever- 
sal(T) . Now applying the parity transformation eq.(6) to the time reversed 
S— matrix we get 





(PT)S(PT)- 1 



"Demanding the invariance of the S— matrix under PT— transformation" we 
get Sll = Srr implying = Ír. It may also be remarked that in the case 
when the Hamiltonian is both parity and time reversal invariant (separately) 
then P invariance alone implies Ír = t L and r R = r L , while T invariance alone 
leads totR — tL- This may be compared to what is obtained in the case where 
we have violation of P and T with PT conservation. 

The reflection amplitudes, however, for left and right moving partides are 
different and are given by 

[q 2 - (k + ~\) 2 ]Sin(qa) e~ ika 

r L = i = . (8) 

2qk Cos(qa) — i (q 2 + k 2 — X 2 ) Sin(qa) 

[q 2 — (k — \) 2 ]Sin(qa) e~ ika 

tr = i • (9) 

2kq Cos(qa) — i (q 2 + k 2 — X 2 ) Sin(qa) 
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Fig. 1. The figures above are for Model-I. The sòlid lines are for left incident 
beam and the dotted ones are for right incidence. l(a): Variation of binding energy 

(/3 = y^fPo with the strength of the imaginary part (À = ^pr) of the potential. 

Here we have taken 2 ^j / ° = 1 and a = 1. l(b): Variation of transmission coefficient 
(|i| 2 ) with the wavenumber (k) for left (or right) incident beam. l(c): Variation of 
reflection coefficient (|r| 2 ) with the wavenumber (k) for left (or right) incident beam. 
l(d): Variation of the deviation from unitarity (|r| 2 + \t\ 2 — 1) with the wavenumber 
(k) for left (or right) incident beam. 

For concreteness consider vàlues of the potential parameters to be range a = 
10, real part of potential given by = 100 and the strength of the imaginary 
part by = 5. We depict the energy variation of the transmission coefficient 
(which is the same for left and right incident beams) as a function of the wave 
number k in Fig.l(b). For the same potential the reflection coefficients for 
right and left incident beams are depicted in Fig.l(c). The quantity |r| 2 + |í| 2 
which is equal to unity for real potentials will now depart from that value 
because of the imaginary part and this is a measure of the extent of inelasticity. 
Accordingly we plot |r| 2 + \t\ 2 — 1 for both left and right incident beams in 
Fig.l(d). 



As a second example of a P-Pseudo Hermitian Hamiltonian we choose a system 
of a particle in one dimension governed by 

H=£- v 5(x) + i\[5(x - a -) - S(x + a -)}. (10) 
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Fig-2(c) Fig.2(d) 

Fig. 2. The figures above are for Model-II. The sòlid lines are for left incident 
beam and the dotted ones are for right incidence. 2(a): Variation of binding energy 

(/3 = \J^j]jr-) with the strength of the imaginary part (À = ^pr) of the potential. 
Here we have taken ^pa = 1 and range a = 1. 2(b): Variation of transmission 
coefficient (|í| 2 ) with the wavenumber (k) for left (or right) incident beam. 2(c): 
Variation of reflection coefficient (|r| 2 ) with the wavenumber (k) for left (or right) 
incident beam. 2(d): Variation of the deviation from unitarity (|r| 2 + |í| 2 — 1) with 
the wavenumber (k) for left (or right) incident beam. 

Here again = PHP~ l but unlike case I the imaginary part of the potential 
is not proportional to the derivative of the real part. Defining 2i p il = fi, 
^pr — A and 2j |^ = f3 2 where B is the binding energy, the eigenvalue condition 
obtained from the jump conditions for the wavefunction becomes 

8/3 3 - 4fi(3 2 = À 2 (l - e^ a )[fi{l - e~ pa ) - 2/3(1 + e pa )]. (11) 



Disregarding the unphysical double root of this equation at j3 — this is solved 
numericaily and the variation of the binding energy is shown as a function of 
the strength of the imaginary part of the potential in Fig. (2a) where we have 
for illustration chosen the real part to be given by fi, — 2i p a = 1 and the 
range a — 1. It may again be noted that there exists no bound state (in the 
usual meaning of the word) below a certain critical value (/i cr ) of the real part. 
This may easily be found from the eigenvalue condition [eg.(9)] by putting its 
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non-trivial root equal to zero and this leads to 

4À 2 a ,_ N 
4 + \ 2 a 2 



The transmission and reflection amplitudes for model-II are given by 



R 



i-*Í) + (è) 2 [-(i 



qe ika + (1 + i£) e 



2ika] 



(13) 



'2k [ 



(i-t + £) + (i-è)[(i + 



; A \ \_ p ika 
' 2k > 2k e 



— C.C. 



l-4) + (è) 2 h(l 



+ (1 + ij-)e 2ífca 



(14) 



'2fc' 



(l + t + #)-(l + è)[(l + 



L 2khk° 



ika 



— C.C.I 



(l-^) + (è) 



^ïtçika 
k 



+ (1 + i|-)e 2ífca 



(15) 



where c.c. stands for the complex conjugate of the preceeding term in the 
braces. 

The reflection and transmission coefhcients for left and right incident beams 
for this model and the departure from elàstic unitarity are shown through 
Fig.(2b) to (2d). For both the models considered it can be easily seen that 
the condition for poles of the reflection and transmission amplitudes is exactly 
the equation for the binding energies, and thus the analytic structure in the 
energy (k 2 ) plane of these amplitudes shows bound state poles. They also have 
a cut along the real axis corresponding to a branch point at k 2 = but the 
discontinuity across the cut is provided by inelastic unitarity arising from the 
imaginary part of the potential. 

Thus we have discussed pseudo-Hermitian Hamiltonians in the context of one 
dimensional complex optical potentials and have shown that the energy eigen- 
values are real for the strength of the imaginary part less than a critical value 
above which we obtain a complex conjugate pair of eigenvalues which create 
some difnculties in their interpretation. We also discuss transmission and re- 
flection from such complex barriers and show that the reflection coeflicient 
for left incident partides is different from that due to those coming from the 
right. We believe it is important to gain experience working with such pseudo- 
Hermitian Hamiltonians as their full physical significance is not yet very clear. 

We gratefully acknowledge help from discussions with Swarnali Bandopadhyay 
and Abhishek Choudhury. 
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